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ABSTRACT 


Tani's  integral  method  is  extended  to  treat  laminar  two-dimensional 
compressible  boundary  layers  with  heat  transfer  and  arbitrary  pressure 
gradient  for  both  attached  and  separated  flows.  A  carefully  chosen 
one-parameter  family  for  the  velocity  profiles  and  a  "universal" 
stagnation  enthalpy  profile  are  assumed  for  attached  flows.  The  . 
accuracy  of  the  method  is  examined  by  comparing  the  results  with 
several  "exact"  numerical  solutions  and  satisfactory  agreement  is 
obtained.  For  separated  flows  one-parameter  families  are  assumed 
for  both  the  velocity  and  stagnation  enthalpy  profiles.  In  this  case 
the  accuracy  of  the  method  is  poor;  however,  suggestions  are  made 
as  to  how  it  might  be  improved  within  the  present  framework. 
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I.  INTRODUCTION 

The  phenomenon  of  flow  separation  is  present  in  almost  all 
fluid  mechanical  devices.  Its  presence  is  seldom  welcomed  since  it 
can  cause  reduced  efficiencies,  increased  drag,  buffeting,  control 
surface  "buzz",  center  of  pressure  shift  on  cylindrical-flare  bodies  , 
and  many  other  troublesome  effects.  In  fact,  in  many  cases  the  ons 
of  separation  puts  an  upper  limit  on  the  performance  --  the  well  known 
"stalling"  of  an  airfoil  is  a  good  example  of  such  a  limit.  On  the  other 
hand,  there  are  situations  where  separated  flows  may  be  beneficial, 
such  as  for  the  reduction  of  drag  and  heat  transfer  at  hypersonic  speeds  . 
Despite  the  obvious  importance  of  flow  separation,  the  problem  has 
escaped  analytic  treatment  because  of  its  complex  nature,  and  remains 
a  poorly  understood  and  essentially  unsolved  problem. 

Separation  may  be  defined  by  introducing  the  concept  of  a 
"limiting  streamline".  Because  of  the  no  slip  condition,  one  cannot 
strictly  speak  of  a  streamline  "at  the  wall".  At  an  infinitesimal  distance 
away  from  the  wall,  however,  the  flow  has  some  finite  velocity  and  hence 
some  definable  direction.  Thus  the  limiting  streamline  is  given  by  the 
limiting  flow  direction  as  the  wall  is  approached.  Since  a  streamline 
cannot  end  in  the  fluid  it  must  either  pass  on  downstream  to  infinity 
or  close  in  the  fluid.  Separation  is  defined  as  the  position  at  which  the 
limiting  streamline  leaves  the  wall  and  enters  the  interior  of  the  fluid. 
Reattachment  is  defined  as  the  position  where  the  streamline  joins 
either  the  surface  or  another  fluid  streamline.  In  two  dimensional 
flow,  the  slope  of  the  limiting  streamline  at  separation  and  at  re- 
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attachment  is  defined  as  (dy/dx)  _q  =  lim  (v/u)  . 

The  point  at  which  the  surface  shear  stress  vanishes  also 
coincides  with  the  separation  point  in  two-dimensional  flow.  For 
three-dimensional  flows,  however,  a  vanishing  shear  stress  is  not  a 
sufficient  condition  for  separation;  thus  the  definition  of  separation 
based  upon  the  limiting  streamline  concept  is  preferable  to  that  based 
upon  the  notion  of  zero  shear  stress. 

*In  the  usual  first  order  boundary  layer  theory  the  required 
"input"  pressure  distribution  is  given  by  the  inviscid  external  flow. 

If  flow  separation  is  present,  however,  the  flow  pattern  and  hence  the 
pressure  distributions  can  be  drastically  altered  from  what  they  would  be 
were  the  fluid  inviscid.  The  present  boundary  layer  methods,  then, 
can  predict  when  separation  is  likely  to  occur;  but  if  separation  does 
occur,  in  many  cases  they  give  little  reliable  information  about  the 

flow  near  the  separation  point  and  in  particular  say  nothing  about  the 

0 

details  of  the  flow  behind  se(>aration.  When  the  flow  is  supersonic,  the 
pressure  field  impressed  by  the  external  flow  is  related  to  the  local 
inclination  of  the  external  flow,  which  in  turn  depends  on  the  "upwash" 
induced  by  the  growth  of  the  boundary  layer.  Thus  the  "feedback  loop" 
is  closed  and  in  this  respect  the  problem  is  somewhat  simpler  than  in 
the  subsonic  case. 

Consider  briefly  the  physical  flow  situation  in  a  typical  super- 

3 

sonic  separating  and  reattaching  flow  .  Sketch  A  on  page  3  shows  the 
flow  and  the  pressure  distribution  in  a  compression  corner.  Typically 
the  separated  flow  region  is  characterized  by  a  more  or  less  constant 
pressure  aft  of  separation  followed  by  a  rising  pressure  just  before 
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SKETCH  A 


reattachment.  At  the  separation  point  the  velocity  along  the  limiting  or 
dividing  streamline  is  zero.  Because  of  mixing,  as  the  fluid  proceeds 
downstream,  the  velocity  along  the  dividing  streamline  increases.  The 
fluid  is  thus  "prepared"  for  the  reattachment  process.  It  is  clear  that 
because  of  the  conservation  of  mass  flow  inside  the  dividing  streamline 
there  must  be  regions  of  reverse  flow.  In  order  for  the  flow  to  reattach 
the  fluid  along  the  dividing  streamline  must  be  brought  to  rest  and  hence 
the  flow  must  experience  a  pressure  rise  prior  to  reattachment.  The 
whole  flow  process  is  a  complicated  interaction  between  the  external 
flow  and  the  viscous  flow  field  —  the  external  flow  adjusts  itself  so  as 
to  affect  che  viscous  region  in  such  a  way  as  to  achieve  reattachment. 

Steady  flows  similar  to  that  just  discussed  also  occur  in  other 
shock  wave  boundary  layer  interactions^,  behind  blunt  based  bodies  in 
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supersonic  flow  ,  and  near  the  leading  edge  of  sharp  nosed  airfoils  . 

The  flows  can  be  either  wholly  lanriinar,  wholly  turbulent  or  "transitional". 

© 

where  transition  takes  place  between  separation  and  reattachment. 

Several  analyses  have  been  devised  to  deal  with  certain  super¬ 
sonic  separated  flows.  The  base  pressure  problem  has  received  a 

great  deal  of  attention  and  approximate  methods  have  been  developed 

3  6 

by  Chapman  .  et  al.  to  treat  the  laminar  case  and  by  Korst  and  his 

co-workers  for  the  turbulent  case.  In  both  methods,  the  details  of  the 

flow  in  the  recirculating  region  are  in  effect  neglected.  The  mixing 

process  is  assumed  to  take  place  at  constant  pressure  and  to  be  the 

same  as  that  which  occurs  in  the  classical  free  boundary  mixing.  The 

analyses  are  valid  when  the  thickness  of  the  boundary  layer  at 

separation  is  zero.  For  these  cases  they  are  found  to  predict  results 

which  compare  favorably  with  experiment.  However,  for  flows  such  as 

the  shock-wave  boundary  layer  interaction  the  upstream  boundary  layer 

is  of  a  size  comparable  to  the  maximum  height  from  the  wall  to  the 

dividing  streamline  after  separation,  and  the  analyses  break  down. 

In  such  cases  the  previous  "history"  of  the  boundary  layer  becomes 

important  and  the  sizable  reverse  flow  velocities  cause  the  mixing 

process  to  depart  from  the  classical  free  boundary  mixing. 

The  usual  one-parameter  Karman- Pohlhausen  momentum 

integral  method  for  attached  boundary  layer  flows  and  its  extension  by 

Thwaites^,  Rott  and  Crabtree®  Cohen  and  Reshotko^  and  others^® 

is  inadequate  for  separated  and  reattaching  flows.  As  shown  in  Sketch 

A  there  are  regions  between  separation  and  reattachment  where  the 

static  pressure  is  very  nearly  constant  and  reversed  flow  occurs  near 
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the  surface.  Since  the  velocity  profile  is  determined  solely  by  the 

local  pressure  gradient  in  the  Karman- Pohlhausen  method,  a 

Blasius-type  profile  would  be  obtained  for  the  pressure  ’’plateau  ’ 

region  if  the  Pohlhausen  quartic  is  employed.  The  reversed-flow 

profiles  found  by  Stewartson^^  along  the  ’’lower  branch”  of  the  Falkner- 

13 

Skan  solution  were  incorporated  into  the  Thwaites  method  by  Curie  , 
but  it  is  not  clear  that  this  special  family  provides  the  required 
flexibility.  Curie's  computed  pressure  distributions  do  not  show  the 
inflection  in  pressure  as  evidenced  by  experiments. 

In  an  attempt  to  ’’unhook”  the  velocity  profile  from  the  local 
pressure  gradient  Crocco  and  Lees^^  introduced  a  new  momentum 
integral  method  in  which  the  profile  is  determined  by  a  single  inde¬ 
pendent  shape  parameter  not  explicitly  related  to  the  pressure  gxadient. 

.  ,  17,  18 

The  Crocco- Lees  method  has  been  developed  quite  extensively  , 

and  has  been  used  to  calculate  such  problems  as  shock-wave  boundary 
layer  interactions,  yielding  fairly  good  results.  The  main  drawback 
of  the  method  is  that  it  relies  upon  a  mixing  coefficient  that  is  not 

well  known  for  separated  flows. 

Some  of  the  arbitrary  features  of  the  Crocco-Lees  method  can 
be  eliminated  by  adopting  and  extending  the  two- moment  method 
developed  by  Tani^^  for  attached  boundary  layers.  In  this  scheme  the 
velocity  profile  is  still  determined  by  a  single  parameter,  say  the  non- 
dimensional  slope  at  the  surface,  a(x),  but  this  parameter  is  independent 

of  the  Thwaites -Pohlhausen  pressure  gradient  parameter  A  (x)  =  - - . 

The  development  of  the  boundary  layer  is  determined  by  integrating  the 
two  simultaneous  first  order  differential  equations  for  a(x)  and  A  (x) 
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obtained  by  taking  the  integrals  of  momentum  and  mechanical  energy 
across  the  layer. 

Lees  and  Reeves^  applied  this  promising  method  to  adiabatic 

separated  flows.  The  purpose  of  the  present  study  is  to  extend  this 

scheme  to  the  more  general  case  of  arbitrary  heat  transfer  at  the 

surface.  Theoretical  analyses  of  heat  transfer  in  separated  regions 

have  been  limited  both  in  scope  and  in  number.  Chapman  examined 

the  effects  of  heat  transfer  and  mass  injection;  however,  the  laminar 

^  ,  .22 

layer  was  required  to  be  thin  and  at  constant  pressure.  Carlson  s 
theory  is  one  of  the  few  that  considered  reverse  flow  velocity  profiles, 
but  this  analysis  too  was  for  constant  pressure.  Curie  treated  the 
shock-wave  boundary  layer  interaction;  however,  when  heat  transfer 
was  present  this  method  gave  rather  poor  results  for  the  prediction  of 
separation, and  its  accuracy  in  the  separated  flow  region  has  not  been 
established. 

It  is  well  known  that  surface  heat  transfer  can  have  a  large 

effect  on  the  behavior  of  attached  boundary  layers.  For  instance,  in  a 

positive  pressure  gradient  cooling  the  surface  delays  separation  and 

heating  the  surface  moves  separation  upstream.  One  of  the  objectives 

of  the  present  investigation  is  to  determine  the  extent  to  which  the 

effect  of  surface  cooling  persists  in  separated  and  reattaching  flow 

regions.  This  effect  could  be  significant  in  determining  the  Mach 

number  along  the  dividing  streamline  and  hence,  the  extent  of  the 

23 

separated- reattaching  flow  itself  . 

In  the  interest  of  simplicity  only  steady  two  dimensional  laminar 
boundary  layers  are  considered  here.  The  flow  up  to  the  separation 


7 


point  is  treated  first  and  results  from  the  approximate  analysis  are 
compared  with  some  "exact"  solutions.  Tani's^^  method  is  generalized 
by  representing  the  velocity  profiles  as  a  weighted  mean  of  the  Blasius 
and  the  average  of  the  non-adiabatic  similarity  separation  profiles, 
rather  than  a  quartic.  The  weighting  parameter  is  again  the  non- 
dimensional  slope  at  the  surface.  The  analysis  is  developed  to  deal 
with  the  flow  beyond  separation  and  one  example  is  briefly  considered. 

Only  cases  for  which  the  external  velocity  is  prescribed  are 
computed.  The  present  study  is,  however,  a  preparation  for  treating 
the  interaction  between  the  viscous  layer  and  the  external  flow. 
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II.  BOUNDARY  LAYER  INTEGRAL  EQUATIONS 

II.  1.  The  Stewartson  Transformation 

The  equations  of  the  steady  laminar  two-dimensional  com¬ 
pressible  boundary  layer  for  perfect  fluids  are: 

Continuity: 


Momentum: 


0  = 


Energy: 


It  is  assumed  that  these  equations  are  valid  for  the  flow  beyond 

separation  as  well  as  for  the  attached  flow.  There  has  been  some 

question  raised  as  to  whether  or  not  the  full  Navier  Stokes  equations 

24 

are  required  at  the  separation  point.  Oswatitsch  demonstrated  that  a 
regular  solution  of  the  Navier-Stokes  equations  exists  in  the  neighborhood 
of  the  separation  point.  *  In  fact,  results  identical  to  his  are  obtained 
if  only  the  usual  boundary  layer  terms  are  kept.  Thus  by  using  integral 
methods  it  should  be  possible  to  pass  through  this  region  without  too 
much  difficulty.  The  assumption  of  negligible  normal  pressure  gradients 
for  the  separated  flow  seems  to  be  a  reasonable  one,  except  possibly  in 


*  The  special  singular  solution  found  by  Goldstein^  may  not  be 
the  one  that  occurs  in  nature. 


certain  cases  near  reattachment  where  the  curvature  of  the  external 
streamlines  is  large. 

For  simplicity  the  viscosity  law  is  taken  to  be 


(5) 


As  shown  by  Cohen  and  Reshotko^,  £qs.  (1),  (2),  and  (3)  may  be  trans¬ 
formed  into  the  form  of  the  two  dimensional  incompressible  laminar 
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boundary  layer  equations  by  means  of  Stewartson's  transformation  , 


even  when  the  flow  is  not  adiabatic. 

A  stream  function  is  defined 


% 


- 


(6) 


and  the  following  variables  are  introduced: 


</y=  ^  #</y 


(7) 

(8) 


The  transformed  incompressible  co-ordinates  are  denoted  by 
upper  case  letters  X  and  Y.  The  subscript  e  refers  to  conditions  at  the 
edge  of  the  boundary  layer,  where  the  flow  is  assumed  to  be  isentropic, 
and  the  subscript  oo  refers  to  conditions  in  the  free  stream.  By  using 
Eq.  (5)  and  the  assumptions  that  Cp  is  constant  and  that  Pr  =  1,  the 
following  equations  are  obtained  for  the  flow  in  the  incompressible  plane; 
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f  =  Ue  ^  (itS)  i  ^ 

^  ay  dx  dY' 


where  S  is  a  dimensionless  enthalpy  defined  by 


(12) 


and  h^  is  the  local  stagnation  enthalpy. 

The  stream  function  has  been  replaced  by  the  transformed 
velocities  (U,  V)  defined  by 


U-  Vr 
V’-K 


(13) 


and  the  resulting  relation  between  the  transformed  and  the  physical 


longitudinal  velocities  is 


u 


CL 

<7# 


£)qs.  (9)  -  (11)  are  subject  to  the  following  boundary  conditions: 


U(X.  0)  =  0 
V(X,  0)  =  0 

S(X,  0)  =  =  constant  (14) 

lim  S  =  0 

Y  -►(» 

lim  U  =  U  (X) 

Y  -»•(» 
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II.  2.  Integral  Form  of  Equations 

When  Pr  =  1  it  seems  reasonable  to  assume  that  the  velocity 
(momentum)  and  the  thermal  boundary  layers  are  of  the  same  thickness. 
Then  integrating  Eqs.  (10)  and  (11)  across  the  boundary  layer  between 
Y  =  0  and  Y  =  -A  and  making  use  of  Eq.  (9)  the  following  momentum  and 


energy  integral  equations  are  obtained 

Ueddi’-  *Z/2*H*i 

5Gc  '  3x 

(15) 

«  dX 

-  Z  4^^  X 

(16) 

Following  Tani^^,  the  momentum  equation,  Eq.  (10),  is  multi¬ 
plied  by  U  and  integrated  across  the  boundary  layer  to  obtain  the  first 
moment  of  the  momentum  equation; 

A 

U.d0i*\  z{3i.z^  =  *3^^“ [(&]%!¥  (‘^1 

dx  »i*  y  dx  lu  y  ' 

0 

The  boundary  layer  characteristics  in  the  incompressible  plane, 
such  as  the  displacement  thickness  6>,  the  momentum  thickness  0^  , 
the  energy  thickness  0.*  ,  the  enthalpy  thickness  ^  ,  and  the  enthalpy 
flux  thickness  ^  ,  are  defined  as  follows: 

o 

o 


(19) 


(20) 


27 

Eqs.  (15)  -  (22)  have  been  given  in  similar  form  by  Foots  . 
Through  Stewartson's  transformation  the  various  boundary 
layer  characteristic  thicknesses  in  the  physical  plane  may  be  related 

a 

to  those  in  the  incompressible  plane.  Thus 

i  *  mtSiJ  (23) 

%'Mr  -  «•*  (25) 

**  <» 


where 
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III.  VELOCITY  AND  TOTAL  TEMPERATURE  PROFILES 

Following  Tani^^,  Poots^^  expressed  the  velocity  and  total  tem¬ 
perature  profiles  as  fourth  degree  polynomials.  For  each  profile, 
four  of  the  coefficients  were  determined  by  fitting  the  boundary  con¬ 
ditions  at  the  wall  and  the  edge  of  the  boundary  layer.  The  remaining 
coefficient  (identified  with  the  gradient  at  the  wall)  for  each  layer  was 
used  to  characterize  the  shape  of  the  profile.  Thus  the  velocity  and 
temperature  profiles  are  each  members  of  a  separate  one-parameter 

family.  The  development  of  these  profiles  along  the  surface  and  the 
o}  dU 

growth  of  0^  (or  — i—  )  is  found  by  integrating  the  three 

V  00 

simultaneous  differential  equations,  Eqs.  (15)  -  (17)  . 

In  the  present  paper  it  is  also  assumed  that  the  profiles  can  be 
expressed  as  members  of  one-parameter  families.  However,  the  pro¬ 
file  shapes  were  determined  in  a  somewhat  different  manner,  because 
the  use  of  Tani's  quartic  for  the  velocity  profile  was  found  to  lead  to 
large  errors  near  separation  for  the  case  of  the  cold  walls.  (This 

point  is  discussed  further  in  Section  IV.  1.  1.  ) 

28 

Cohen  and  Reshotko  present  similar  Falkner-Skan  type 
solutions  (i.  e. ,  when  =  CX^)  for  the  laminar  compressible  boundary 
layer  with  heat  transfer.  When  the  separation  profiles  in  the  trans-  , 
formed  plane  for  various  values  of  are  normalized  and  compared, 
these  velocity  profiles  do  not  collapse  to  one  universal  curve  (Figure  1). 
Thus  it  would  appear  that  at  least  two  parameters  are  required  to 
represent  the  profiles  for  the  general  case  of  an  arbitrary  .  The 
addition  of  a  second  parameter  would  require  the  addition  of  another 
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differential  equation  (for  example  the  second  moment  of  the  momentum 
equation)  in  order  to  solve  the  flow  problem.  *  However,  to  avoid  the 
added  complexity  of  a  second  parameter,  a  kind  of  "mean"  one- 
parameter  velocity  profile  to  be  used  for  all  was  chosen  in  the 
following  way: 

The  velocity  profile  is  written  as  =  f(Y/4  )  +  a(X)  g(YA  )  » 

where  a  is  identified  with  the  gradient  at  the  wall.  Thus  when  a  =  0, 
U/Ug  =  f(YA)  ,  u  e.  ,  the  separation  profile.  The  function  f(Y/A  )  was 
determined  by  taking,  in  a  sense,  the  "average"  of  the  exact  similarity 
separation  profiles  from  Reference  28  for  various  values  of  S^.  The 
representative  average  profile  chosen  for  f(Y/^  )  was  the  same  as  the 
exact  separation  profile  for  =  -0.  8,  except  for  slight  modification 
in  order  that  the  boundary  layer  have  finite  thickness. 

For  a  =  a„-  ,  the  velocity  profiles  in  the  trans- 

formed  plane  for  all  values  of  reduce  to  the  Blasius  profile.  Thus, 
a  g(YA)  =  (U/U  -f(Y/!i)  ,  where  the  subscript  BL  refers  to 

the  "Blasius"  values.  Since  "a"  corresponds  to  the  gradient  at  the  wall, 
the  function  g(Y/A)  was  then  determined.  Explicitly, 


/  8(U/UJ„l\ 

'  'yA=o 

By  definition  f'(0)  =  0,  and  g  was  chosen  such  that  g'(0)  =  1.  For  the 
profiles  as  chosen,  a^j^  =  1.99.  The  functions  f(Y/^  )  and  g(Y/^)  are 


*  Based  on  the  similar  solutions,  this  second  parameter  could 
instead  be  determined  as  a  function  of  and  thus  the  additional 
differential  equation  would  not  be  required.  This  approach  would  make 
the  tabulation  of  the  boundary  layer  functions  defined  by  Eqs.  (31)  -  (35) 
very  involved. 


® 
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given  in  Table  1  and  are  shown  in  Figure  2.  This  one -parameter 
family  is  used  to  describe  the  velocity  profiles  in  the  separated  flow 

region  as  well  as  for  attached  flow. 

28 

Cohen  and  Reshotko  present  total  temperature  profiles  in 
the  incompressible  transformed  plane  for  various  values  of  and 
pressure  gradient  parameter  p.  If  these  profiles  are  "normalized"  by 
"scaling"  the  normal  distance  from  the  wall  so  that  all  profiles  have  the 
same  gradient  at  the  wall,  upon  comparison  an  interesting  result  is 
obtained.  (Figure  3  shows  the  profiles  at  separation  for  various 
compared  with  the  "flat-plate"  profile.  )  As  long  as  the  flow  is  attached 
the  S/S^  profiles  can  be  represented  quite  well  by  one  "universal"  curve. 
This  universal  curve  is  given  by  Crocco's  integral  of  the  energy  equation 
for  the  flat  plate,  namely,  S/S^  =  1  -  (U/Ug)^^,  where  (U/Ug)gj^ 

is  the  "Blasius"  profile.  Thus,  for  attached  flow  the  thermal  profile 
is  taken  as  S/S^  =  “  {^(Y/a  )  +  glY/A^^  •  In  the  separated 

flow  region  the  thermal  profiles  can  no  longer  be  represented  by  this 
"universal"  profile.  The  separated  thermal  profiles  are  assumed  to  be 
the  one- parameter  family, 

S/Sw  =  -  MY/a  )^  +  b(X)  g(Y/A  ) 

where  for  convenience  f(Y/A)  and  g(Y/A  )  are  the  same  functions  as 
those  used  for  the  velocity  profiles. 

Summarizing,  the  velocity  profile  for  both  attached  flow  and 
separated  flow  (as  long  as  the  height  of  the  reverse  flow  region  is  not 
too  large)  is  taken  as 

U/Ug  =  f(Y/A)  +  a(X)  g(Y/A) 

For  attached  flow  the  "universal"  thermal  profile  is  used 


(28) 
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S/S^  -  1.  99  8(Y/a)  ,  (29) 

i.  e. ,  b  =  -1.  99,  while  for  separated  flow  the  thermal  profile  is  given  by 

S/S^  =  j^-f(YA)^+  b(X)  g(Y/A)  (30) 

Now  that  the  profile  shapes  are  decided  upon,  Eqs.  (18)  -  (22) 
can  be  integrated  graphically  to  give  the  various  non-dimensionalized 
boundary  layer  thicknesses  in  terms  of  a  and  b.  Thus 

(6.*)M  =  .4204  -.0651a  =  D  (31a) 

(0.)/Z^  =  .09080  +  .02616a  -  .00842a^  =  E  (31b) 

0.*/A  =  .  1368  +  .0360a  -  .00655a^  -  .001 182a^  =  F  (31c) 

/A 

( VUg^)  / (  aU/aY)^  dY  =  1.  763  -  .  5040a  +  .  2068a^  =  (Q/4F)  (32a) 


o 

=  .4204  +. 0651b  =  W  .  E  (32b) 

^/S  A  =  .  0908  +  .0456a +.  01947b +.  00842ab  (32c) 

=  J  .  E  =  Z  .  F 

where  J  =  (^/S^O.)  and  Z  =  (£/S^0*)  ; 

H  =  D/E  ,  G  =  F/E  (33) 

(20./Ug)  (au/aY)Y^Q  =  P  =  2aE  (34) 

-(lH/sJ')  (as/aY)Y_o  =  R  =  -2b  J  .  E  =  -2bz  .  F  (35) 


The  quantities  given  by  Eqs.  (31)  -  (35)  are  functions  of  a  and  b 
only  and  their  numerical  values  are  tabulated  in  Table  2. 

Eqs.  (15)  -  (17)  are  now  rewritten  in  the  form 

Vt  -f  =v^P 

y.  .  zc.'9:^(3*  zs„%)4p!  .  ^  q 

dX  ^  (37) 


(36) 
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dX  ^  (38) 

Given  the  external  velocity  wall  temperature  and 

the  initial  conditions,  the  set  of  first  order  differential  equations, 

Eqs.  (36)  -  (38),  can  now  be  solved  for  the  three  unknown  functions 
a(X),  b(X),  and  0.(X)  . 


e 
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IV.  SOLUTIONS  OF  THE  BOUNDARY  LAYER  INTEGRAL  EQUATIONS 
IV.  1.  Attached  Flow 

For  attached  flow,  the  "universal”  S/S^  profile  is  used  and  thus 
"b"  is  constant  and  numerically  equal  to  -1.  99.  The  problehu  is 
simplified  and  reduced  to  solving  the  two  first  order  differential 
equations,  Eqs.  (36)  and  (37). 

IV.  1.  1.  Similarity  Solutions 

When  the  flow  is  a  similar  Falkner-Skan  type  flow  (i.  e. ,  when 
Ug  =  CX^)  a  =  const,  and  Eqs.  (36)  and  (37)  reduce  to  two  algebraic 
equations.  Because  of  the  way  in  which  the  velocity  and  temperature 
profiles  were  chosen  the  errors  in  the  solution  of  the  integral  equations 
are  largest  for  a  similar  flow  which  is  always  on  the  verge  of  separation 
(i.  e. ,  a  =  0).  For  this  case  values  of  the  pressure  gradient  parameter, 

(3  =  (2m/m+l),  were  calculated  for  various  values  of  S^.  Figure  3 

shows  a  plot  of  p  .  vs.  calculated  by  the  present  method 

™  28 
compared  with  exact  solutions  obtained  by  Cohen  and  Reshotko  . 

It  should  be  mentioned  that  initially  Tani's  quartic  was  chosen 

for  the  velocity  profile  and  a  one-parameter  cubic  was  chosen  for  the 

thermal  profile.  The  full  set  of  equations,  Eqs.  (36),  (37),  and  (38), 

were  solved  for  Pg^paration  ®w  '  "" 

Figure  4,  and  it  can  be  seen  that  when  these  profile  shapes  are  used 

the  integral  method  is  in  considerable  error  for  cold  walls.  This 
error  is  mainly  caused  by  the  large  difference  between  Tani's  quartic 
velocity  profile  and  the  exact  profiles  at  separation.  For  this  reason 
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the  profile  shapes  were  determined  as  stated  in  Section  III,  resulting 

in  improved  accuracy  for  the  cold  wall.  The  qualitative  shape  of  the 

exact  fi  .  vs-S  curve  is  not  matched  by  the  integral  method. 

Kgeparation  w 

Near  =  -1.  0  the  exact  curve  has  negative  curvature  whereas  with 
only  a  one-parameter  family  for  the  velocity  profiles  the  integral  method 
gives  positive  curvature  everywhere.  However,  the  present  method 
gives  values  fairly  close  to  the  exact  solution  except  at  S  =  -1.  0. 

The  displacement,  momentum  and  enthalpy  thicknesses  and  the 
gradient  of  S/S^  at  the  wall  were  calculated  by  the  present  method  for 
the  separation  profile.  These  data  compared  with  the  exact  results  from 
Reference  28  are  shown  in  Table  3. 

The  comparison  is  favorable  except  for  the  displacement  thick¬ 
ness  at  =-1.0  (highly  cooled  wall)  where  the  errors  brought  about 
by  the  one-parameter  velocity  profile  show  up  rather  strongly. 

Again  it  is  repeated  that  the  errors  in  the  present  method  will 
most  likely  be  greatest  for  this  case  of  "incipient  separation". 

IV.  1.  2.  Flow  with  Linearly  Decreasing  External  Velocity  in 
Transformed  Plane 

Calculations  have  been  carried  out  for  the  case  of  a  linearly 
decreasing  velocity  in  the  incompressible  plane,  i.  e.  , 

where  is  the  velocity  at  X  =  0  and  L  is  some  characteristic  length. 

A  Pohlhausen  type  parameter  is  introduced 


(40) 
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Using  Bis.  (39)  and  (40),  after  some  manipulation  Eqs.  (36)  and 
(37)  may  be  put  in  the  following  form  suitable  for  numerical  integration 


1 

II 

[iA- 0] 

(41) 

_  / 

W  0-t) 

if  '<*7 

(42) 

where 

=  x/l 

(43) 

a 

=  (EG^a  -  ^) 

(44) 

=  (1  +  2  Z  -  H 

-  S  W  ) 
w  ' 

(45) 

Jr 

=  (6  +  4  Z) 

(46) 

Now  o,  and  J*  are  functions  only  of  a,  Eqs,  (41)  and  (42)  may 
be  solved  by  eliminating  X  and  numerically  integrating  the  single  result¬ 
ing  differential  equation  in  the  A  ~  tn  Q  plane,  and  then  by  a  simple 
quadrature  transforming  to  the  X  plane,  A  solution  may  also  be 
obtained  by  numerically  solving  Eqs.  (41)  and  (42)  simultaneously  and 
this  second  method  was  used  here. 

Eqs.  (41)  and  (42)  are  subject  to  the  initial  condition  that 

yl  =  0  at  X  =  0  (47) 

The  numerical  integration  was  started  by  the  Runge-Kutta 

29  19  27 

method  and  continued  by  Milnes  method  .  As  Tani  and  Foots  found 
with  a  uniformly  retarded  external  velocity  it  was  difficult  to  carry  out 
the  solution  right  up  to  the  separation  point  because  of  the  rapid  growth 
of  (Hfrir./djK)  near  separation.  (The  behavior  near  separation  is  discussed 
in  Section  V.  1.  )  However,  the  solution  was  carried  out  to  X  =  0.  56  and 
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extrapolated  from  there  to  separation. 

From  the  numerical  computation,  values  of  a  and  A  are 

obtained  for  various  values  of  X.  Using  this  information  the  boundary 

layer  characteristics  in  the  incompressible  plane  such  as  6^*  ,  8^  ,  8^*  , 

<ir  .  X.  (9S/8Y)y^0  (0U/8Y)y^q  can  be  calculated. 

Z7 

Using  the  Hartree-Womer sley  method  Foots  obtained  an  "exact" 

numerical  solution  for  the  above  case  of  S  =1.0  and  a  linearly  decreasing 

external  velocity  (taking  Uj  =  1  and  L  =  8  to  simplify  the  numerical 

calculations).  Foots  also  presented  the  results  of  an  integral  method 

which  amounted  to  solving  the  set  of  three  differential  equations, 

Eqs.  (36)  -  (38)  ,  using  Tani's  quartic  for  the  velocity  profile  and  a 

similar  quartic  for  the  total  temperature  profile.  Calculations  by  the 

,  .  r  «  .  27 

present  method  compare  favorably  with  the  exact  solution  of  Foots  as 
shown  in  Figures  5  and  6.  The  present  integral  method  is  somewhat  less 
accurate,  but  also  simpler  than  the  integral  method  of  Foots.  It  is  also 
expected  that  the  present  method  might  be  more  accurate  than  the 
integral  method  of  Foots  for  the  more  interesting  case  of  a  cold  wall 
since  the  present  velocity  profiles  are  probably  more  realistic  than 
those  represented  by  Tani*s  quartic. 

For  comparison  purposes,  the  boundary  layer  properties  for  the 
case  of  a  cold  wall  with  the  wall  temperature  equal  to  the  initial  tem¬ 
perature  of  the  external  stream,  i.  e.  ,  =  -0.  762  were  computed  by 

the  present  method.  These  results  are  shown  in  Figures  7  and  8. 

The  present  method  for  a  non-adiabatic  wall  predicts  that  the 

13 

heat  transfer  rate  at  separation  is  finite.  Analyses  such  as  Curie's 
which  express  the  temperature  profiles  as  power  series  of  the  velocity 
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ratio,  U/Ug  ,  incorrectly  yield  zero  heat  transfer  at  separation  and 
can  also  reverse  the  sense  of  the  wall  temperature  gradient  in  the 
reverse  flow  region. 

IV.  2.  Separated  Flow 

For  separated  flow  it  is  no  longer  possible  to  use  the  "universal  ' 
S/S^  profile  and  the  full  set  of  differential  equations,  Eqs.  (36),  (37),  and 
(38)  must  be  solved. 

IV.  2.  1.  Similarity  Solutions 

Again  for  the  similar  type  flows,  Eqs.  (36),  (37),  and  (38)  reduce 
to  a  set  of  simultaneous  algebraic  equations.  These  equations  were 
solved  for  the  particular  case  of  =  -0.  8  and  j3  =  -0.  10.  The  velocity 
and  temperature  profile  were  computed  and  are  compared  in  Figure  9 
with  the  exact  solution  of  Cohen  and  Reshotko  .  The  various  integral 
thicknesses  were  not  calculated  in  Reference  28  and  thus  no  comparison 
is  made.  However,  as  can  be  seen  from  Figure  9  the  comparison  of 
the  profile  shapes  is  rather  poor  and  the  present  integral  method  gives 
only  very  rough  estimates  of  such  things  as  skin  friction  and  wall  heat 
transfer  rates.  The  reason  for  these  discrepancies  lies  in  the  inability 
of  the  assumed  form  of  the  velocity  profile  to  match  the  exact  profile 


shape. 
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V.  DISCUSSION  AND  FUTURE  WORK 


The  present  method  points  out  rather  simply  some  interesting 
features  that  occur  near  the  separation  point.  Eliminating  (dO^  /dX) 
from  £qs«  (36)  and  (37)  one  obtains 


Near  separation,  as  a  ^  0  ,  G(a)  goes  through  a  minimum, 
i.  e. ,  (dc/da)  =  0.  Thus  when  (da/dX)  is  finite  and  (do/da)  =  0  ,  a 
unique  value  is  obtained  for  ^  ,  i.  e.  ,  .A  =  A  •  However,  for 
example,  in  the  case  of  a  uniformly  retarded  external  velocity  (Section 
IV.  1.2.)  when  Eqs.  (41)  and  (42)  are  integrated,  it  is  found  that  a 
value  of  A  is  reached  before  separation  such  that  A  <  A  (algebraically). 
From  Eq.  (41)  it  can  be  seen  that  (dA/dX)  <  0  for  all  X  up  to  separation. 
Thus,  when  (dC/da)  =  0  ,  A  A  °  which  implies  that  (da/dX)  =  -«  , 
and  a  singularity  occurs  at  this  point.  It  is  found  that  when  (da/dX)  at 
separation  is  infinite,  (d6*/dX)  is  also  infinite;  however,  (dO./dX)  and 
(de^*/dX)  remain  finite.  Examining  the  wall  shear  stress 


and  since 


-r-a 


y/)r»0 


Ve 

A 


Uc  +  Q:  dUe  _  l^eO.  dA 

7^  dx  A  3x  3x 


(49) 


^  1 

dx  o- 
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[42^]  da. 

dx  (50) 

Keeping  first  order  terms  near  separation,  it  can  be  shown 
easily  from  Eq.  (48)  that 

^  <^separation  '  (8“/5y)y,o  sep^vs^tion  " 

as  assumed  by  Goldstein^^  near  the  separation  point.  As  Goldstein 

found,  for  this  special  type  of  external  velocity  distribution  (when 

(d^Ug/dX^)  =  0  )  the  solution  cannot  be  continued  downstream  of  the 

separation  point. 

30  31  32 

Prandtl  and  later  Meksyn  *  showed  that  the  pressure 

distribution  in  the  region  of  separation  cannot  be  chosen  arbitrarily 

but  must  satisfy  certain  conditions  compatible  with  the  reverse  flow 

region  downstream  of  separation.  Prior  to  separation  A  must  go 

2  /  2  „ 

through  a  minimum  and  near  separation  (d  U^/dX  )  >  0  .  This 
condition  is  evident  from  the  integral  form  of  the  equations,  if  one 
takes  Eq.  (40)  and  examines  the  conditions  for  (dA/dX)  to  change  sign, 
in  order  that  A  pass  through  /i  =  J\,  °  at  separation. 

While  the  special  class  of  flows  where  (d^U^/dX^)  =  0 
(right  up  to  the  separation  point)  lead  to  a  singularity  at  separation  and 
cannot  be  carried  downstream,  exact  solutions  for  such  cases  are 
nevertheless  of  interest  for  checking  approximate  methods. 

@ 
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V.  2.  Velocity  Profiles  in  Separated  Region 

Examination  of  the  solutions  of  the  Falkner-Skan  equation  shows 
that  along  the  lower  branch,  the  maximum  backflow  velocity  is  zero  at 
separation,  reaches  a  maximum  as  p  increases  algebraically  and  then 
decreases  to  zero  again  as  p  ■  ^  0.  The  displacement  thickness 
increases  without  limit  as  p  — ►O  along  the  lower  branch  (i.  e.  ,  at 
p  =  0  the  profile  is  the  classical  free  boundary  mixing  problem  for  zero 
pressure  gradient).  Clearly  such  a  behavior  cannot  be  reproduced  by 
the  kind  of  one- parameter  velocity  profile  chosen  in  the  present  paper. 
The  need  for  such  a  behavior  is  illustrated  by  some  work  of  Reeves* 
for  the  shock-wave  boundary  layer  interaction  on  an  adiabatic  flat  plate. 
Reeves  used  Tani's  quartic  for  the  velocity  profile  and  included  a 
third  integral  moment  equation  to  relate  the  pressure  gradient  to  the 
displacement  effect  of  the  boundary  layer.  It  was  found  that  the 
pressure  did  not  level  off  into  the  usual  "plateau"  region  but  reached 
a  maximum  and  then  decreased  before  rising  again  at  reattachment. 

The  velocity  profile  was  such  that  the  displacement  thickness  could 
not  grow  fast  enough  to  obtain  the  pressure  "plateau". 

Thus  it  would  appear  that  for  the  reverse  flow  region  it  may  be 
necessary  either  to  use  a  two- parameter  velocity  profile  or  to  use 
two  or  more  layers.  Another  simpler  and  promising  method  has  been 
suggested  by  Professor  L.  Lees.  A  one- parameter  family  of  velocity 
profiles  could  be  constructed  based  upon  the  "lower  branch"  Falkner- 
Skan  solutions.  It  should  be  noted  that  it  is  not  necessary  to  relate  the 


♦  private  communication 


parameter  to  any  pertinent  physical  quantity,  but  only  that  the  profiles 
be  denoted  by  this  single  parameter  in  such  a  way  that  the  integral 
properties  may  be  tabulated  (as  in  Table  2), 

y,  3,  Interaction  Between  Viscous  Flow  and  External  Stream 

For  all  the  examples  computed  here  the  external  velocity 
gradient  was  assumed  to  be  given.  In  a  problem  such  as  the  shock 
•wave  boundary  layer  interaction  the  external  velocity  is  not  known 
a  priori  and  the  interaction  between  the  viscous  flow  and  the  external 
stream  must  be  determined.  The  following  equation  is  obtained  by 
applying  the  Stewartson  transformation  to  the  continuity  equation  and 
integrating  across  the  boundary  layer 


where 


<s> 


streamline  direction  angle  relative  to  a  flat  wall 
(oriented  in  the  free  stream  direction)  at  y  =  6, 


e  tan  & 


m 


00 


U« 


And  for  example,  when 


<^  <  <  I,  tan  (9>  c:  C5>  in  Eq.  (51)  and  the 


linearized  Prandtl- Meyer  equation  gives 


where 

M^  =  M  +  £•  :  jef  <<  lA 

e  00  ’ll  00 

Thus  Mg  takes  the  form  of  a  dependent  variable  when  Eqs.  (51) 
and  (52)  are  added  to  the  set  of  Eqs.  (36)  -  (38)  . 
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VI.  CONCLUDING  REMARKS 

Tani's  two-moment  integral  method  has  been  extended  to  treat 
non-adiabatic  two-dimensional  compressible  boundary  layers.  The 
assumption  of  a  "universal"  stagnation  enthalpy  profile  for  all  pressure 
gradients  and  wall  temperatures  is  found  to  be  quite  accurate  for 
attached  boundary  layers  and  provides  a  useful  simplification.  The 
accuracy  of  the  integral  nriethod  is  found  to  be  sensitive  to  the  choice  of 
the  velocity  profile.  By  use  of  the  universal  temperature  profile  and  a 
carefully  chosen  one-parameter  velocity  profile  the  problem  is  reduced 
to  solving  two  first  order  ordinary  differential  equations  when  the 
pressure  gradient  is  prescribed.  Predictions  of  the  boundary  layer 
properties  and  the  separation  point  by  this  method  compare  favorably 
with  "exact"  numerical  solutions. 

'  Flow  beyond  the  separation  point  is  briefly  considered.  The 

"universal"  temperature  profile  is  no  longer  applicable..  With  the 

assumption  of^one-parameter  families  for  temperature  and  velocity 

profiles,  it  is  necessary  to  solve  three  first  order  ordinary  differential 

equations.  By  comparing  the  present  results  with  the  reverse-flow 

28 

Falkner-Skan  profiles  found  by  Cohen  and  Reshotko  one  concludes 
that  the  separated  flow  velocity  profiles  in  any  integral  method  must 
be  described  either  by  a  two-parameter  family,  or  by  the  Falkner- 
Skan  family  itself.  Another  possibility  is  to  use  a  multi-layer  method. 
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FIGURE  1 

SIMILARITY  SEPARATION  VELOCITY  PROFILES 
IN  THE  INCOMPRESSIBLE  PLANE 
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FIGURE  3 

SIMILARITY  TOTAL  TEMPERATURE  PROFILES 
IN  THE  INCOMPRESSIBLE  PLANE 
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THERMAL  BOUNDARY  LAYER  CHARACTERISTICS 
IN  INCOMPRESSIBLE  PLANE  FOR  S„.  =  -  0.  762 


